Complex evolution of the electronic structure from polycrystalline to monocrystalline 

graphene: generation of a new Dirac point 
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First calculations, employed to address the properties of polycrystalline graphene, indicate that 
the electronic structure of tilt grain boundaries in this system display a rather complex evolution 
towards graphene bulk, as the tilt angle decreases, with the generation of a new Dirac point at the 
Fermi level, and an anisotropic Dirac cone of low energy excitations. Moreover, the usual Dirac point 
at the K point falls below the Fermi level, and rises towards it as the tilt angle decreases. Further, 
our calculations indicate that the grain-boundary formation energy behaves non-monotonically with 
the tilt angle, due to a change in the the spatial distribution and relative contributions of the bond- 
stretching and bond-bending deformations associated with the formation of the defect. 
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Graphene - an isolated layer of sp^-bonded carbon 
atoms arranged in a honeycomb structure - was until very 
recently, a "theoretical" reference system for the study 
of the properties or the "real" sp^-bonded carbon forms, 
such as fuUerenes, nanotubes, and graphite 
Since the recent report of the isolation of a stable single- 
atom-thick carbon layer, by exfoliation of graphite [5|, 
graphene itself has occupied the center stage of Mate- 
rials Physics, as a paradigmatic system for "relativis- 
tic" condensed-matter phenomena, as well as a promising 
material for Nanoelectronics, due to its exceptional elec- 
tronic properties. Graphene is a null-gap semiconductor 
with a vanishing density of states at the Fermi level, and 
electronic bands that are linear and isotropic within 
eV from the Fermi level. This linearity and the presence 
of two sub-lattices imply that charge carriers in graphene 
effectively behave as massless chiral "relativistic" parti- 
cles, being described by Dirac's equation. Due to 
the chiral nature of the electronic excitations, which leads 
to the absence of backscattering, graphene holds ballistic 
charge transport on the microscale, even at room temper- 
ature, and with high concentrations of defects and impu- 
rities. P, 0, Q In the last few years, scientific interest 
in graphene has rapidly intensified, and the material is 
expected to play a major role in Nanoelectronics in the 
future. 

Presently, common synthesis routes for graphene 
are the original exfoliation method, Q that produces 
monocrystalline graphene samples, and epitaxy, mostly 
on SiC [7] and metallic substrates Q. Recent works on 
epitaxially-grown graphene report the occurrence of su- 
perstructures interpreted as Moire patterns 0, based 
on STM, AFM, and STS measurements. Moire pat- 
terns and superstructures, associated with layer stack- 
ing, as well as the occurrence of bulk and surface grain 
boundaries, are topics of prominence in the physics of 




highly-oriented pyrolytic graphite (HOPG) itself. [11 
12 . 13 . 14, [l^ Grain boundaries (GB) are among the 



FIG. 1: (Color online) Grain boundary geometry: the relative 
orientation between grains is defined by the angle between 
vectors Ml and Mr. The distance between the two grain 
boundaries in the supercell is d, and L is the period of the 
pent agon- heptagon pattern along the grain boundary. 



most commonly occurring extended defects in HOPG, be- 
cause of its polycrystalline character [l3| . For large scale 
graphene production and application, it is expected that 
synthesis methods will be epitaxy based, being quite con- 
ceivable that polycrystalline samples will be produced. 
Indeed, the occurrence of GBs on the graphene layer 
has been recently reported, speculated as a probable 
source of l ong -range electronic perturbations in graphene 
on Si02, [lO| and further, the lower carrier mobility of 
epitaxial graphene (when compared with exfoliated sam- 
ples), in a macroscopic-size graphene field-effect transis- 
tor, has been tentatively assigned to the electronic per- 
turbations associated with GBs. [3] 

Previous works have addressed the electronic proper- 
ties of disordered graphene, but these have focused on 
the effects of point defects, corrugation, and extended 
edge states in graphene ribbons, drawing only specula- 
tive conclusions in what regards the electronic states of 
GBs. fl*] In this scenario, addressing properties of GBs in 
graphene is of primary relevance. In the present work, we 
employ first principles calculations to examine the ener- 
getics and electronic properties of GBs in graphene. We 
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TABLE L Geometric parameters and formation energy per 
unit length Ef (in eV/A), of GB supercells containing A'^^ 
atoms, a is the tilt angle and L is the GB period (in A), as 
indicated in Fig. [T] E^*^ and i^ef"'* are the contributions to 
the Keating-model total elastic energy per unit length E"^' 
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focus on the structural GB model proposed by Simonis et 

al. 



Ill ] who observed a large-angle tilt boundary on the 
surface of HOPG, on STM experiments. These authors 
proposed that, in the absence of stress, the observed GB 
consists of a periodic structure that can be described as a 
regular succession of pentagon-heptagon pairs, as shown 
in Fig. [1] Based on this model, we investigate three GBs 
of different periodicities along the boundary, hence with 
different relative orientation between the grains. 

We find that, while large-angle tilt GBs do not in- 
troduce localized states at the Fermi level in graphene, 
various resonance peaks appear in the density of states 
of the material, over a broad energy range, starting at 
energies of ~0.3 eV from the Fermi level, in agreement 
with recent experimental work. [l5| More importantly, 
the changes in electronic structure with the GB tilt angle 
indicate a non-trivial evolution towards graphene bulk, as 
the GB angle decreases: we observe the generation of a 
new Dirac point at the Fermi level, which lies on a line 
that evolves towards the T-M direction of the graphene 
Brillouin zone (BZ), with the usual Dirac point at the 
K point falling below the Fermi level, and rising towards 
it as the tilt angle decreases. The Dirac cone of low- 
energy excitations around this new Dirac point is non- 
isotropic, with the effective "speed-of-light" depending 
on the direction in the BZ away from the Dirac point. 
Furthermore, our calculations indicate that, within the 
structural model we consider, the GB formation energy 
does not behave monotonically with the period of the GB 
(or, equivalently, with the tilt angle). A Keating anal- 
ysis of the elastic energy, associated with the formation 
of the defect, indicates that a change in the the spatial 
distribution and relative contributions of the stretching 
and bending deformations leads to the non-monotonicity 
indicated by our results. 

All calculations are performed using Kohn-Sham den- 
sity functional theory Il7|, the generalized-gradient ap- 
proximation (GGA) [l8| for the exchange-correlation 
term, and norm-conserving Troullier-Martins pseudopo- 
tentials [l^, to describe the electron-ion interaction. 
We use the LCAO method implemented in the SIESTA 
code [i^l , with a double-zeta pseudo-atomic basis set plus 
polarization orbitals, with an energy cutoff of 0.01 Ry. 
Structural optimization is performed until the total force 
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FIG. 2: Band structure and density of states for (a) GBl, (b) 
GB2, and (c) GB3 geometries. The Fermi level is indicated 
by the dashed line. Brillouin-zone lines are shown in Fig. [3] 



on each atom is less than 0.02 eV/A. In order to simulate 
an isolated honeycomb sheet, we use supercells that are 
periodic along the graphene plane, and are surrounded by 
a 33 A vacuum region, such that the interactions between 
each layer and its periodic images are negligible. 

Periodicity along the graphene plane requires the su- 
percell to contain two GBs of opposite tilt angles (a GB 
and the corresponding "anti-GB"), as shown in Fig.[T_ 
the GBl geometry. The experimental value in Ref. 
for the relative orientation between grains is 21°, defined 
here by the angle a between vectors and Mr, drawn 
respectively on the left and right grains adjacent to the 
GB, as shown in Fig. [T] This model can be extended to 
GBs with smaller tilt angles, by adding lines of hexagons, 
such that the period L of the pentagon-heptagon pattern 
along the GB increases. In the L — > oo (a — > 0) limit, we 
recover the perfect single-crystal graphene lattice. We 
study three different GB geometries, with the theoretical 
values for L and a indicated in Table H] GBl is the model 
proposed in Ref. [11], with a = 21.8° and L = 6.6 A; GB2 
has a = 13.3° and L = 10.9 A; and for GB3 a = 9.6° 
and L — 15.2 A. In order to ensure that we simulate the 
properties of an isolated GB, we consider supercells with 
increasing distances d between the GBs and their peri- 
odic images. Formation energy results are converged for 
d = 14.9 A. The geometric parameters a, d, and L are 
shown in FiglU with values for a and L given in Table [H 

The band-structure of the GBl {a = 21.8°), is shown 
in Fig. [2ja). The Brillouin zone for this GB supcrccU 
is shown in Fig. ^h). The F— Y direction is along the 
GB, and the F— X direction is perpendicular to it. In 
Fig. [3I^b), we also indicate lines we denoted as F— M and 
F— K, which correspond to these high-symmetry direc- 
tions in a single-crystal graphene sheet with the orienta- 
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TABLE II: Fermi velocities (in lO^m/s) for the GB Dirac 
electron cones, along the indicated directions. 



FIG. 3: (Color online) Electronic energy dispersion for poly- 
crystalline graphene. (a) The anisotropic Dirac cone around 
the D point, for the GBl. (b) The Brillouin-zone direc- 
tions and regions for the GBl. (c) and (d) The color-coded 
constant-energy isolines (with the energy increasing from vi- 
olet to red) for the GBl and GB3, respectively, (e) and (f) 
The wedge-shaped dispersion around the K-X line shown in 
(b). 



tion of the grain on the right side of the GB. F— M (and 
F— Y) and F— K (and F— X) converge to the correspond- 
ing directions in the a — > limit, where a single-crystal 
graphene sheet is recovered. A mirror plane, perpendic- 
ular to the graphene sheet and through the geometric 
center of the GB, relates the F— M and F— K lines of the 
two adjacent grains. Note that the Fermi level or Dirac 
point occurs on the F— Y line, in the point marked as D 
in Fig.[3]Jb), and that band crossings, which lie below the 
Fermi level, occur at the X and K points. 

Thus, we observe the generation of a new Dirac 
cone 
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on the F— Y direction, which is anisotropic, 
with a Fermi velocity (the effective "velocity of light") 
that depends on the direction from the Dirac point in 
k-space. This Dirac cone is shown in Fig. EJa), for the 
region around the D point indicated in Fig. [Sljb)). The 
energy isolines are shown Fig. [3fc) , where the anisotropy 
of the cone is clearly displayed: the more (less) elongated 
is the isoline along a given direction, the smaller (larger) 
is the Fermi velocity along that direction. The band- 
structures for the GB2 (a = 13.3°) and GB3 (a = 9.6°) 
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are shown in Figs.[2tb) and (c), respectively. Note again 
the occurrence of a Dirac point at the Fermi level, on 
the F— Y line, in both cases. Both display anisotropic 
Dirac cones as well. The isolines for the GB3, shown 
in Fig. EJd), are less anisotropic than in the GBl case. 
In Table HI] we include the Fermi velocities for the Dirac 
cones of the three GBs, along the indicated directions. 

The nature of the electronic dispersion in the region 
around the X and K points, also indicated in Fig. ^h), 
are shown in Fig. ^e) and (f), for the GBl and GB3, 
respectively. For the GB3, there is very little dispersion 
along the X-K direction, resulting in a straight wedge 
cutting through this line, while in the GBl, the shape of 
the energy dispersion is also wedge-like, but a "flowery" 
shape develops at higher energies. Furthermore, the band 
crossings at the X and K points move up towards the 
Fermi level, as a decreases, which is consistent with the 
a limit, where the Fermi level occurs at the K point. 
The energy difference between the Dirac point at D and 
the K point is 1.00 eV for the GBl, 0.73 eV for the GB2, 
and 0.64 eV for the GB3. 

These results indicate a very complex evolution of 
the band-structure of polycrystalline graphene with the 
GB angle. Since fivefold and sevenfold topological de- 
fects constitute the core of low-energy dislocations in 
graphene, [S^] they are probably ubiquitous in any re- 
alistic model of GBs in this material. Hence, the oc- 
currence of a new Dirac point along the GB direction 
may prove a robust feature of polycrystalline graphene. 
We note that previous works have found that, while 
pentagon-heptagon pairs and the related SW defect in- 
troduce no resonant states at the Dirac point in graphene, 
resonance peaks appear in the density of states starting 
at a few tenths of an eV from the Fermi level, a result 
that we have also reproduced with the ab initio method 
employed here. Our GB electronic structure calculations 
show, however, that the periodic superstructure formed 
by these dislocations along the GB line shares with per- 
fect graphene the vanishing gap and the Dirac-like na- 
ture of electronic excitations, but in a rather complex 
structure, with direction-dependent Fermi velocities, the 
generation of a new Dirac point on the GB direction, and 
a wedge-like dispersion around the K point that should 
evolve towards the graphene Dirac cone as a ^ 0. We 
expect charge transport in polycrystalline graphene to 
reflect the anisotropic structure of the Dirac cones we 
obtain in our calculations. 
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FIG. 4: (Color online) Spatial distribution of bond-bending 
and bond-stretching energies of grain boundaries in graphene. 
(a) Bend for GBl. (b) Stretch for GBl. (c) Bend for GB3. 
(d) Stretch for GB3. 



We turn now to the GB energetics. The GB for- 
mation energy per unit length is given by Ef — 
(Egb - Ebuik) f^L, where Egb and Ebidk are the cal- 
culated total energies of the GB and bulk graphene su- 
percells, respectively. The factor of two on the right-hand 
side accounts for the presence of two GBs in the cell. 

Results in Table |T] show that Ef varies non- 
monotonically with the GB period. Among the three 
GBs, GBl with a tilt angle of 21.8° has the smallest for- 
mation energy, while GB2 with a = 13.3° has the highest 
Ef value. We expect that i?/ — > as a ^ 0. Our re- 
sults indicate that Ef initially increases, before reaching 
this asymptotic limit. Given the absence of broken bonds 
in the GB geometries, Ef is primarily of elastic nature. 
In order to understand this non-monotonic behavior, we 
use a Keating model to analyze the bond-bending i?gf"'^, 
and bond-stretching E^f^ contributions to the GB elas- 
tic energy E^i ■ The results are shown in Table HI For 
the GBl we observe that E^f^"^ is a factor of 3.2 larger 
than E^f^, while for the GB2 and GB3 structures the 
j^bend I j^strch pg^yg dccreascs to 2.1 and 1.8 respectively. 
These results remain essentially unaltered for different 
choices of the ratio between the stretching and bending 
energy parameters of the Keating model. Moreover, the 
spatial distribution of Sgf""* and E^^^'^'^^ changes as a 
decreases, as shown in Fig. El For the GBl, -E^f"'^ is 
largely concentrated on the atoms at the very core of the 
boundary, mainly in the pentagon-heptagon pairs, while 
i?!*'' is a little more spread out towards the interior of 



the grains. For GB3, we see that both £;°f"'^ and Ell"" 
have significant contributions from atoms in the interior 
of the grain, mostly on the hexagon lines that at the 
pentagon-heptagon at the boundary. Thus, by concen- 
trating the elastic energy on bending distortions at the 
defect core, the GBl structure relaxes to a lower energy 
state than GBs with smaller angles, which leads to the 
non-monotonic behavior oi Ef. 

To conclude, ah initio calculations indicate that the 
electronic structure of tilt GBs in graphene display a 
complex evolution towards graphene bulk, as the GB an- 
gle decreases, with the generation of a new Dirac point at 
the Fermi level, lying on a line that evolves towards the 
F-M direction of the graphene Brillouin zone, at the ver- 
tex of an anisotropic electronic-energy cone. Moreover, 
the usual Dirac point at the K point falls below the Fermi 
level, and rises towards it as the tilt angle decreases. Fur- 
thermore, our calculations indicate that the GB forma- 
tion energy behaves non-monotonically with the tilt an- 
gle, due to a change in the the spatial distribution and 
relative contributions of the bond-stretching and bond- 
bending deformations associated with the formation of 
the defect. 
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